Thermodynamic transitions in inhomogeneous rf-wave superconductors 
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We study the spectral and thermodynamic properties of inhomogeneous d-wave superconductors 
within a model where the inhomogeneity originates from atomic scale pair disorder. This assumption 
has been shown to be consistent with the small charge and large gap modulations observed by 
scanning tunnelling spectroscopy (STS) on Bi2Sr2CaCu2 08+cc- Here we calculate the specific heat 
within the same model, and show that it gives a semi- quantitative description of the transition 
width in this material. This model therefore provides a consistent picture of both surface sensitive 
spectroscopy and bulk thermodynamic properties. 
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Introduction. There is accumulating evidence from 
STS that at least some of the families of superconduct- 
ing cuprate materials are "intrinsically" inhomogeneous 
at the nanoscaleP, 0, H, 0, in the sense that all sam- 
ples of the given material exhibit electronic disorder over 
length scales of ~ 25A which cannot be removed by 
any standard annealing procedure. In particular, differ- 
ential conductance maps of Bi2Sr2CaCu208+x(BSCCO) 
reveal the existence of spectral peaks reminiscent of the 
coherence peaks of a homogeneous d-wave superconduc- 
tor whose energy varies by a factor of 2-3 over the sam- 
ple. There is currently widespread interest in this phe- 
nomenon, but no consensus as to its origin, nor even as 
to whether the gaps between these peaks may be taken 
to correspond directly to the local superconducting order 
parameter in the sample, or whether some are related to 
a second, competing order 0. While these experiments 
have been proposed as evidence of a general nanoscale 
inhomogeneity in cuprates, there have been several criti- 
cisms. One class of objections points out that other ma- 
terials, particularly YBa2Cu307_5(YBCO) exhibit much 
narrower NMR linewidths, suggesting that in this mate- 
rial at least, the electronic inhomogeneity cannot be a 
pronounced property of the bulkH- A stronger critique 
has been offered by Loram et ai.Q even with regard to 
the BSCCO material. These authors point out that if 
one associates with each "gap patch" a local doping cor- 
responding to the hole density required to produce a gap 
of that size in the average phase diagram, then the dis- 
tribution of the doping would be of order Ap ~0.1. Such 
a large range of local doping levels would in turn pro- 
duce a large inhomogeneous distribution of regions with 
different Tc's, sufficient to yield a transition width ATc 
of order the average Tc itself. Since this is contrary to 
specific heat measurements |^ 0, these authors conclude 
that the inhomogeneity is either a surface phenomenon, 
or represents a distribution of antinodal scattering life- 
times rather than of order parameters. Clearly it is very 
important to determine whether results from surface sen- 
sitive probes reflect the bulk properties of these materials. 



Recently, localized resonances were imaged by STS at a 
bias of -960 meV and identified as the O interstitials pri- 
marily responsible for doping the BSCCO mater ialfioj. 
This work also noted the strong positive correlation be- 
tween the positions of these dopants and the local spec- 
tral gap, and furthermore argued that the charge inho- 
mogeneity in the sample was considerably smaller than 
anticipated on the basis of previous scanning tunnelling 
topo grap hs integrated to smaller bias voltages. Nunner 
et ai.jTl| then argued that these and a number of other 
experimental observations could be explained by the 
simple hypothesis that the disorder caused by the dopant 
atoms was primarily in the Cooper rather than density 
channel, i.e. that each dopant modulated the BCS pair 
interaction on an atomic scale. 

In this Letter we argue that the lack of strong charge 
modulations within the pair disorder model of Ref. |ll| 
also allows one to avoid the argument of Loram et al. Q 
regarding the transition width. The model, with parame- 
ters chosen to reproduce the gap maps and other correla- 
tions of the STS data on BSCCO, is shown to yield results 
within mean field theory which are consistent with a rel- 
ative specific heat transition width, ATc/Tc^ of ~ 20%, 
despite the fact that at T = gap modulations of order 
100% are observed. The observed widths must there- 
fore be attributable to pair disorder alone since the sam- 
ples are homogeneous at the mesoscale (as determined, 
e.g., by optical and scanning electron microscopy). Pre- 
dictions of the theory for the formation and persistence 
of superconducting islands near the transition should be 
verifiable by high-temperature STS measurements. 

Model. There are many studies of inhomogeneous 
super conductivity 0, only few of which are directly rel- 
evant to the questions addressed below. If one adds non- 
magnetic disorder to an ordinary s-wave superconduc- 
tor, one expects only small changes in superconducting 
one-particle properties such as the gap and Tc due to 
Anderson's argument [l^. Nevertheless, as disorder in- 
creases and the mean free path i becomes comparable to 
the Fermi wavelength A^, localization effects can destroy 
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superconductivity. Recently, local aspects of this transi- 
tion were studied by Ghosal et al. jl5j , who showed using 
numerical solutions of the Bogoliubov-de Gennes (BdG) 
equations including nonmagnetic disorder that near the 
transition the highly disordered system separates into is- 
lands with finite order parameter surrounded by an in- 
sulating sea. In the d-wave case, where ordinary dis- 
order is pair breaking and the transition occurs when 
~ ^0, where is the coherence length, several numer- 
ical studies have inve stig ated local properties of disor- 
dered systems[ll[ll[li,[l3. Only Ref. [H considered 
the finite temperature transition, however, pointing out 
that the inhomogeneity induced by ordinary disorder (in 
this case the suppression of the order parameter around 
strong scattering centers) strongly affects the transition, 
leading to a Tc onset (defined by the formation of su- 
perconducting islands as temperature T is reduced from 
above), considerably higher than the Tc predicted by the 
usual disorder- averaged theory of Tc suppression, anal- 
ogous to the theory of Abrikosov and Gorkov[2(j. The 
width of the transition was not discussed explicitly in this 
work, nor was the effect of pairing disorder considered. 

More recently, it has been proposed that Tc can even 
be increased by local inhomogeneities compared to its 
value for the homogeneous system [i^, 113 • Within a 
weak coupling BCS framework it was shown that peri- 
odic modulations of the pairing and/or electron density 
can lead to a substantial enhancement of Tc when the 
characteristic modulation length scale is of order ^oH^- 
This agrees with recent numerical calculations for the at- 
tractive Hubbard model with various disorder distribu- 
tions of the interaction, verifying that Tc enhancement by 
inhomogeneity is a general result arising from the prox- 
imity effect 123. These calculations did not investigate 
the details of the transition below the Tc where the local 
order parameter A first nucleates. 

In the following model calculation we use the standard 
(i-wave BCS Hamiltonian 

H=Y^ ti^-4^c^-^+^(l/i-/i)4^Q^+^ (^Aijcl^c]^ +H.C.) , 

{ij)a icr (ij) 

(1) 

where in the first term we include nearest t and next- 
nearest f = —0.3t neighbor hopping. For the chemical 
potential /i, we set fi = —t in order to model the Fermi 
surface of BSCCO near optimal doping ~ 16%. Xl(ij) 
denotes summation over neighboring lattice sites i and j. 
Disorder is included by the impurity potential Vi = Vofi 
where fi = ^^e'x.-p{—ris/X)/ris^ r^g being the distance 
from a defect s to the lattice site i in the plane. Dis- 
tances are measured in units of >/2a, where a is the Cu- 
Cu distance. Note that the particular Yukawa form of 
fi is merely a convenient way to vary the smoothness of 
the potential landscape through the parameter A. The 
(i-wave order parameter Aij = gij{ci^Cji — Cj|Q|)/2 is 
determined self-consistently by iterations of 

= ^Y1 (^n{i)vn{j) + Vn{i)un{j)) tanh(^), (2) 

n 



until convergence is achieved. Here, {En^Un^Vn} is the 
eigensystem resulting from diagonalization of the BdG 
equations associated with Eq.Q. The pairing interac- 
tion gij varies in space relative to its homogeneous value 
go = 1.16t as gij = go-\-Sg{fi-\-fj)/2^ where Sg is the mod- 
ulation amplitude and i^j are nearest neighbors. In the 
following, when including potential (rs channel, Vq 7^ 0) 
or pair (ri channel, Sg 7^ 0) disorder, we make sure to 
adjust fi and go, respectively, in order to maintain the 
same doping and average gap as the corresponding ho- 
mogeneous system. Note that in this approach the in- 
clusion of spatial pair potential variations is purely phe- 
nomenological. However, a likely candidate for the pair- 
ing modulations is the oxygen dopants. Various possible 
microscopic origins for this unusual disordered state have 
been discussed in Refs. 0, E3- There it was also argued 
that Ts disorder alone cannot account for the STM data. 

Results. Fig. Q shows, in a field of view (FOV) of 
49nm x 49nm, the experimental LDOS at -960 meV (a), 
and the corresponding gap map (b) obtained by McEl- 
roy et al. . The bright resonances in Q^a) reveal the 
location of the oxygen dopants. In Fig. ^c) we show 
the impurity potential generated by using each of these 
dopants (833 in the above FOV) as a defect located out 
of the Cu02 plane. This map resembles (a) as it should. 
The experimental FOV is modeled as a 90 x 90 lattice 
system rotated 45 degrees with respect to the 3.83 A 
long Cu-Cu bonds, i.e. it includes 2 x 90 x 90 sites and is 
aligned with the experimental FOVj^. The theoretical 
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FIG. 1: (Color online.) (a) Experimental dl/dV map [arb. 
units] at -960 meV of an optimally doped BSCCO sample from 
Ref. [Io|. (b) Experimental gap map [meV] in the same region 
as (a), (c) The theoretical impurity potential extracted from 
(a) assuming distance from plane = 0.5 and A = 0.5. (d) 
The gap map resulting from using (c) as the 'pairing potential' 
in the BdG equations with 6g = 1.5t. Both (c) and (d) are 
shown in units of t. 
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gap map (as extracted from the LDOS) using (c) as a ti 
potential with = 1.5t at T = is shown Fig. dd). We 
find that gap maps consistent with experiment are found 
for 1.0 < 5g/t < 2.0. The correlation coefficient between 
the gap maps^b) and^d) is 0.17 26j, a reasonably large 
value given the simplicity of the pure n calculation and 
the fact that the correlation between the dopant positions 
and the experimental gap map is 0.3|Tnj. 

Now we turn to the discussion of the results for T 7^ 0. 
Below, we set = 50 (corresponding to the upper left 
5/9th square of each image in Fig^ and study the re- 
sulting order parameter (OP) maps, entropy 5'(T), and 
specific heat C{T) with focus on the temperature region 
near Tc. The number of iterations necessary for con- 
verged results increases dramatically near Tc, whereas 
just a few degrees away from Tc we typically find that 
25-50 iterations suffice. In Fig. [2 we show the OP map 
Ai = {Ai^i^x + \,i-x - \,i+y - \,i-y) /4, for temper- 
atures near Tc. The bulk transition temperature for the 
homogeneous system is Tco = 0.182t. Here, one clearly 
sees how separate islands of finite start forming above 
Tco when cooling down and eventually overlap at lower T. 
In principle, STS measurements at T close to Tc should 
be able to probe these superconducting islands. 

In order to address the effect of the inhomogeneity on 
the transition widths, we first calculate the quasiparticle 
entropy S{T) according to the well-known expression 

^ = -2^ [/(^Oln/(^n) + /(-^n)ln(/(-^n))], (3) 

where f{E) is the Fermi distribution function. From 
the resulting entropy S{T) curve we extract the specific 
heat at constant volume C{T) by the usual expression 
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FIG. 3: (Color online) Specific heat C{T) for the impurity 
distribution of Fig. QJc) vs T: (a) clean 50 x 50 BdG system 
(solid, black line), "infinite" clean system (dashed, black line), 
and Ti disordered with 5g — l.Ot (thick solid, blue line); 5g — 
1.5t (thick dashed, red line); 5g = 2.0t (thick dash-dotted, 
black line), (b) Comparison of C{T) for the Nambu channels 
n and T3 with otherwise identical impurity parameters. 



C{T) = T{dS/dT). Fig. Elja) compares the electronic 
specific heat for the clean and pair disordered systems 
described by different disorder strengths Sg. In Fig. [S^b) 
we compare C{T) for the different Nambu channels ri 
and Ts. We can check the BdG results by comparing to 
the pure "infinite" system result (black, dashed line in 
Fig. OJa)) obtained by replacing En 



FIG. 2: (Color online) OP maps, parameters from Fig. 1(d): 
T = 0.18t (a), T = 0.20t (b), T = 0.22t (c), and T = 0.24t (d). 



with ek the normal state dispersion with hopping inte- 
grals and chemical potential identical to the values given 
above, and Ak = Ao{T){coskx — cos ky)/ 2 with Ao(T) 
solved selfconsistently from the usual gap equation. It is 
clear that the 50 x 50 BdG lattice problem is large enough 
to capture the sharp transition width of the clean system. 

Conventional potential disorder corresponding to ~ 
10% dopant atoms causes only negligible broadening in 
the specific heat transition width compared to the pair 
disordered case. This is shown in Fig. Ofb), where both 
curves are produced from the impurity distribution in 
Fig. dc). The transition is sharp because potential scat- 
terers cause local suppressions of the order parameter, 
whereas in interstitial regions A decreases with temper- 
ature in a manner similar to the pure system. We have 
also studied distributions of strong scatterers at the per- 
cent level, and find similarly sharp transitions. On the 
other hand, potential disorder smooth on a scale of can 
lead to large modulations of fi which, even within BCS 
theory, will give larger transition widths. In this case, 
however, the associated local spectra are not consistent 
with STS measurements, as shown in Ref. 0. 

As evident from OJa), pairing disorder with param- 
eters fixed to yield reasonable variations of the gap 
size and graininess at low temperature ^ leads to a 
broadened transition width similar to the experimental 
observations J] . This is shown more clearly in Fig. ^ 
where we compare the experimental specific heat C{T) 
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FIG. 4: (Color online) (a) Specific heat transitions vs T for 
YBCO and BSCCO optimally doped powders, from Ref. 0. 
(b) C(T) — Cn{T) (N: normal state) vs T with parameters 
from FigOIa), assuming that Tco = 0.182t equals Tc = 90K. 



for clean YBCO and optimally doped BSCCO near 
with our results for the clean and pair disordered case. 
It is clear that in the case of YBCO, the transition is ex- 
tremely narrow, and the deviations from the mean- field 
treatment discussed here are consistent with weak 3DXY 
critical fluctuations over a range of a few degrees [2^ 12^. 
On the other hand, in the case of BSCCO, it has been 
clear for some time that the 3DXY description of the spe- 



cific heat transition is not appropriate, despite the fact 
that according to a naive application of the Ginzburg cri- 
terion 3DXY fluctuations should be more visible 8], and 
it has in fact been discussed as closer to a Bose-Einstein 
condensate (BEG) with specific heat exponent of a = — 1, 
although measured values are closer to a = —0.70 • This 
interpretation remains controversial, however. Other au- 
thors have suggested that a 3DXY divergence in C{T) cut 
off by bulk nanoscale inhomogeneities over a small range 
ATc/Tc ^ 0.03 could be consistent with the dataf^ [2£^. 
Here we have put forward a similar scenario, but our mi- 
croscopic model implies that the inhomogeneities domi- 
nate over a larger range ATc/Tc ~ 0.1 — 0.2 to be con- 
sistent with STS. Additional contributions from critical 
fluctuations must be added to the mean field effects dis- 
cussed here to obtain a complete description. 

Conclusions. We have presented theoretical calcu- 
lations for d-wave superconductors with atomic scale 
pair disorder, using impurity parameters appropriate to 
reproduce semi-quantitatively the gap maps produced 
by STM experiments on optimally doped BSCCO, and 
shown that the experimental specific heat transition 
in this system can be explained by this model as well. 
This suggests that substantial nanoscale electronic in- 
homogeneity is characteristic of the bulk BSCCO system. 
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